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ABSTRACT 



Context. While low mass clouds have been relatively well studied, our picture of high-mass star formation remains unclear. Infrared Dark 
Clouds appear to be the long sought population of cold and dense aggregations with the potential of harbouring the earliest stages of massive 
star formation. Up to now there has been no systematic study on the temperature distribution, velocity fields, chemical and physical state toward 
this new cloud population. 

Aims. Knowing these properties is crucial for understanding the presence, absence and the very potential of star formation. The present paper 
\Q aims at addressing these questions. We analyse temperature structures and velocity fields and gain information on their chemical evolution. 
{*""*) ■ Methods. We mapped the (/, A") = (1, 1) and (2,2) inversion transitions of ammonia in 9 infrared dark clouds. Our observations allow the most 
^5 ' reliable determination of gas temperatures in IRDCs to date. 

Results. The gas emission is remarkably coextensive with the extinction seen at infrared wavelengths and with the submillimeter dust emission. 
Our results show that IRDCs are on average cold (T < 20 K) and have variations among the different cores. IRDC cores are in virial equilibrium, 
O are massive (M > 100 M Q ), highly turbulent (1-3 km s -1 ) and exhibit significant velocity structure (variations around 1-2 km s -1 over the 
£j cloud). 

ffi ' Conclusions. We find an increasing trend in temperature from IRDCs with high ammonia column density to high mass protostellar objects and 
. . hot core/Ultracompact Hn regions stages of early warm high-mass star formation. The linewidths of IRDCs are smaller than those observed in 
r*- . high mass protostellar objects and hot core/Ultracompact Hn regions. On basis of this sample, and by comparison of the ammonia gas properties 
within a cloud and between different clouds, we infer that while active star formation is not yet pervasive in most IRDCs, local condensations 
might collapse in the future or have already begun forming stars. 

Key words. infraredTSM: Clouds, ISM: Molecules, ISM: Structure, Molecular Data, Radio Lines: ISM 
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1. Introduction 

Infrared dark clouds (IRDCs) are cold, dense molecular clouds 
seen in silhouette against the bright diffuse mid-infrared (MIR) 
emission of the Galactic plane. They were discovered dur- 
ing mid-infrared i maging surveys wit h the Infrared Space 
Observatory (ISO, IPerault etaD Il99fj|) an d the Mid-course 
Space Experiment (MSX. lEgan et alJl 998). 

In an initial census of a ~180° long st rip of the G alactic 
plane (between 269°< / < 91°, b = ±0.5°1. lEgan et al.ldl998h 
found ~2000 compact objects seen in absorption against bright 
mid-infrared emission from the Galactic plane. Examination of 
2MASS, MSX and IRAS images of these objects reveals that 
they appear as shadows at all these wavelengths, although they 
are best identified in the 8.3/im MSX band, because, first, the 
7.7 and 8.6 /mi PAH features associated with star-forming re- 
gions contribute to a brighter background emission and, sec- 
ond, the MSX 8.3/im band is more sensitive than the satel- 
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lite's other bands. Rece ntly, we have reviewed the obser va- 
tional studies on IRDCs (Me nten. Pillai. & Wvrowskil200^) . 

While low mass clouds have been relatively well studied, 
our picture of hig h-mass star formation remains unclear (see 
lEvans et "all 12002). IRDCs appear to be the long sought pop- 
ulation of cold and dense aggregations with the potential of 
harbouring the earliest stages of massive star formation. It is 
likely that some of the stars forming in them are massive (lu- 
minosities of submm conden sation s range up to 10 4 L Q ). A 
recent study by lOrmel et al.l J2005I) on an IRDC towards the 
W5 1 Giant Molecular Cloud (GMC) suggests that sources of 
~ 300 M Q are embedd ed within the cores, most likely proto- 
stars. Recently, we (Pill ai et all20 05) reported a detailed study 
of the strongest submm peak in the IRDC Gl 1 . 1 1 -0. 12, where 
we find clear evidence of a heavily embedded protostar. 

The salient results on IRDCs are summarised below. The 
IRDCs observed so far have sizes of 1-10 pc and have mostly 
a filamentary morphology. On t he basis of LV G calculations 
of mm H 2 CO observations ICarev et alJ Jl998) find that typi- 



2 



T. Pillai, F. Wyrowski, SJ. Carey and K.M. Menten: Ammonia in Infrared Dark Clouds 



cal IRDCs have gas densities of n > 10 6 cm -3 and temper- 
atures of T < 20 K. Kinematic distances determined from 
H2CO observations using a standard Galactic rotation curve 
that ranges between 2.2 and 4.8 kpc dCarev et alJ 1998) indicate 
that the clouds are not lo cal. All observed IRDCs in the sam- 
ple of lCarev et aD {2000) contain 1-4 bright sub-millimeter 
(submm) dust continuum emission peaks (> 1 Jy/ 14 " beam at 
850 //m) surrounded by an envelope of emission which matches 
the morphology of the IRDC in mid-infrared extinction. The 
cores corresponding to the brightest submm peaks have masses 
100 and 1200 M , except for two clouds in the Cygnus region 
tha t have masses arou nd 40 M . 

iHennebelle et alJ J200 lh in a systematic analysis of the 
ISOGAL images extracted about 450 IR DCs, for which they 
derive 15//m optical depths of 1 to 4. lTevssier et alJ J2002I) re- 
ported that Large Velocity Gradient (LVG) model calculations 
of HC3N, 13 CO, and C ls O yield densities larger than 10 s cm -3 
in the densest parts. The authors claim to find kinetic tempera- 
tures between 8 and 25 K based on CH3CCH observations; the 
higher values being found toward embedded objects, however 
a detailed analysis is hitherto unpublished. 

Observations of molecules in IRDCs so far have concen- 
trated on just a few species. Up to now there has been no sys- 
tematic study on the temperature distribution, velocity fields, 
chemical and physical state toward this new cloud population. 
Knowing these properties is crucial for understanding the pres- 
ence, absence and the very potential of star formation. The 
present paper aims at addressing these questions. We have 
started with a survey of a sample of 9 IRDCs in (1,1) and 
(2,2) cm rotational transition s of ammonia (NH 3). This sam ple 
has been studied before bv lCarev et al . ( 1998) and Carev et al. 
(2000) and was selected on the basis of the large extent and 
high contrast of the IRDCs against the MIR background. 

Ammonia has proven to be an important t ool in measur- 
ing th e physical conditions in molecular clouds fao & TownesI 
1983i|. Since only the lowest NH3 energy levels are expected 
to be populated in cool dark clouds (T < 20 K), their physi- 
cal conditions can be probed using the (1,1) and (2,2) inversion 
transitions in the metastable (J, K) rotational levels of ammo- 
nia. Radiative transitions between different Zf-ladders are for- 
bidden, therefore the lowest levels are populated only via colli- 
sions. The optical depth can be determined from the ratio of the 
hyperfine satellites. Thus, the population of the different lev- 
els can be estimated and hence the temperature of the gas de- 
termined. In addition, recent chemical models reveal that NH3 
(and also N2H + ), does not deplete from the gas phase for the 
densi ties observed in IRDCs (< 10 6 cirT 3 ) (B ergin & La nger 
1997). Thus NH3 is an excellent tracer of the dense gas where 
many other molecules would have heavily depleted. 

In §|2] we give details of our observations with the 
Effelsberg 100m telescope. In {j3] we discuss the data reduc- 
tion and present the correlation between gas emission and MIR 
absorption. In ^4] we derive the rotational temperature, gas ki- 
netic temperature and NH3 column density. Furthermore, we 
analyse the velocity structure, estimate dust mass, virial mass 
and NH3 abundance. In ^5] we compare the core gas properties 
(temperature, NH3 linewidths and column density) with that of 
other populations of objects that are thought to trace the early 



stages of high mass star formation. We do a similar comparison 
with local dark clouds well-studied in NH3. Finally, we spec- 
ulate on a possible formation mechanism of IRDCs involving 
supernova remnants. 

2. Observations 

We mapped the IRDCs listed in Table. ^ with the 
Effelsberg 100m telescope of the Max-Planck-Institut fur 
Radioastronomie in October 1999. The frontend was the fa- 
cility 1.3 cm maser receiver tuned to a frequency of 23.7 GHz 
centered between the NH3 (1,1) and (2,2) transitions. The spec- 
trometer was a 8192 channel auto-correlator used with 2 sub- 
units of 20 MHz bandwidth each. The resulting spectral reso- 
lution was * 0.2 kms _I after smoothing the data to increase 
the signal-to-noise ratio. The beamwidth at the frequencies of 
the NH3 lines is 40" (FWHM). The observations were con- 
ducted in frequency switching mode with a frequency throw 
of 7.5 MHz. The maps toward all sources covered the extinc- 
tion seen in the MSX images and were made with half power 
beamwidth (HPBW) spacing (40"). Alternate scans were inter- 
spaced at half the full beamwidth resulting in * 20 x -\[(2) spac- 
ing. Pointing was checked at roughly hourly intervals by means 
of continuum drift scans on nearby pointing sources. We found 
the pointing to be accurate to within 12". 

Absolute calibration is not an issue in determining the rota- 
tional temperature, because it is solely governed by the ratio of 
the NH3(1,1) hyperfine lines and the (1,1) and (2,2) brightness 
temperatures. However, in order to estimate the excitation tem- 
perature and the column density the data needs to be calibrated. 
The calibration procedure is documented on our webpage 
The important steps are mentioned below. The 100m data (nor- 
mally in CLASS format) is in arbitrary noise tube units (T N ) 
and has to be converted to main beam brightness temperature 
units. We observe a standard flux calibrator with known flux S a 
at wavelength A and thus estimate the conversion factor from 
Tn scale to Jy. 

The main beam brightness temperature Tmb of the calibra- 
tor for a given beamwidth 6 and wave length A is then given 
by 



(7WK) 



(S,,/Jy) x (i/cm) 2 
(0/arcmin) 2 X2.65' 



(1) 



The primary flux calibration is based upon continu um scans 
of NG C7027 assuming a flux density of « 5.6 Jy JOtt et alJ 
11994 . The Ultracompact Hn (UCH11) regions G10.62-0.38 and 
G34.26+0.15 and the quasar J1743-038 were used as sec- 
ondary calibrators which span the entire elevation range and an 
elevation-dependant calibration factor was derived. The max- 
imum rms uncertainty in calibration for elevations lower than 
19° is 15%. 

3. Data Reduction and Analysis 

The spectra were reduced using the CLASS package 
jForveille. Guilloteau. & Lucaslll989l) . For a given source, the 

1 www.mpifr-bonn.mpg.de/staff/tpillai/eff_calib/index.html. 
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spectra were averaged and a polynomial baseline of order 3-5 
subtracted. 

The maps were generated by convolving the original data 
using a Gaussian function with a HPBW of 40". Table^reports 
the coordinates of the central position for each map. 

MSX 8^/m images of these clouds with superimposed con- 
tours of the NH3 (1,1) integrated intensity are shown in Fig.^ 
The velocity range used for integration is within +25 km s~ ! of 
the 'Local Standard of Rest Velocity', Vlsr, hence including 
the satellite lines. The NH 3 maps of G79.27+0.38 and G33.71- 
0.01 have clumpy structures. Some of the clumps may be arti- 
ficial and caused by the low signal-to-noise ratio. 

The reduced (1,1) and (2,2) spectra for the NH3 peak posi- 
tions of the observed sources are shown in Fig. 13 NH3(1,1) 
observations were reduced using "METHOD NH3(1,1)" in 
CLASS to fit the hyperfine structure and derive optical 
depths and linewidths. The standard procedure to analyse the 
NFL (1,1) and (2,2) lines have been described in detail by 
Bac hilleret aD < ll987t) . 

The hyperfine structure of the (2,2) line is too weak to be 
observed; therefore, the optical depth could be obtained only 
for NH 3 (1,1). The NH 3 (1,1) and (2,2) main beam brightness 
temperatures are obtained by fitting the main line with a sin- 
gle Gaussian. The (1,1) and (2,2) linewidth is obtained by the 
hyperfine fitting which accounts for the line broadening due to 
optical depth effects excluding self-absorption. The parameters 
obtained from the fits to the (1,1) hyperfine lines are listed in 
TableEl 

The basic physical parameters, namely the excitation tem- 
perature, rotational temperature, the kinetic temperature and 
ammonia column density, have been derived using the stan- 
dard formulation for NH3 spectra llHo & Townesll983h . Tabled 
summarises the estimates of these parameters toward the cores. 
We give the formal errors (lcr), derived from Gaussian error 
propagation. 



4. Results 

Our sample consists of both 'Infrared dark clouds', with extents 
of ~ 1— 10 pc, and 'Infrared dark cloud complexes', which are 
comprised of multiple individual clouds. IRDC Gil. 11-0. 12 
would thus be a complex (> 4 pc) while G24.63+0.15 would 
be a dark cloud. The IRDCs G79.34+0.33 and G79.27+0.38 
are essentially two parts of the same extended filament in the 
Cygnus-X region connected by a brig ht patch of dust emis sion 
as seen in SCUBA and MSX images (Redm anet al.l2003l) . 



4.1. Cloud Morphology 

There is generally a close match between the ammonia emis- 
sion and the mid-infrared extinction as shown in Fig. ^ 
In G24.72-0.75 and G79.34+0.33, the NH 3 emission peak 
is ho wever correlat e d with bright and compact MIR emis- 
sion. iRedman et alJ J2003I) find that the MIR emission in 
G79.34+0.33 corresponds to a luminous 'Young Stellar Object' 
(YSO) with a strong IR excess which might be interacting with 
the foreground IRDC. No such interaction has been reported 
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Fig. 3. SCUBA 850/im image JCarev et alJl2000h of the cloud 
Gil. 11-0. 12 with NH3 (1,1) integrated intensity as contours. 
The contour levels are multiples of 2<x. SCUBA image has a 
resolution of 14". 

for the IRDC G24.72-0.75 and the nature of the MIR object is 
unknown. 

In general, the cloud geometry is extended and filamentary 
and in no case close to spherical. The mean aspect ratio deter- 
mined by fitting 2D ellipse to the entire NH3 emitting region is 
2.2 and the total extent of the clouds ranges from 0.4 - 8.3 pc. 
Gll.ll-0.12is a filamentary cloud as revealed by the 8/rni ex- 
tinction and the 850//m dust emission as shown in Fig. [3] The 
peaks of the submm emission are strikingly coincident with the 
ammonia cores. The NH3 map reveals at least two peaks to- 
ward the north-east segment and another set of peaks towards 
the southwest extension of the filament suggesting that there 
are several unresolved sub-structures or cores within the ex- 
tended filament. 

The strongest submm emission peak P2 in G28.34 is in the 
close vicinity of the IRAS source 18402-0403. But the peak of 
the ammonia emission for the northern extension seems to be 
offset from that of the IRAS source (see Fig.0}. This could be 
due to the interaction of the IRAS source with the cloud. 

4.2. Source Sizes 

Due to the large distances to the clouds, the compact structures 
within the cloud are marginally resolved with the 40" beam of 
our NH3 observations. Hence, estimating the size of the core 
from NH3 might only deliver upper limits of the order of the 
beam size. We use channel maps to solve the problem of iden- 
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Fig. 1. MSX image of the clouds at 8.3/vm (greyscale) with NH3 (1,1) integrated intensity as contours. The NH3 (1,1) maps 
have a resolution of 40" and the same spacing. The greyscale corresponds to intensity range shown by the wedge in the upper 
right corner (from dark to light). The contour levels are -2cr, lcr, 2<x, 4cr, 6cr... The rms noise levels are given in Table[2] The 
approximate map boundary is also indicated. Tick marks are coordinate offsets (in arcseconds) relative to the positions given in 
Tabled- The bar indicates a length of 1 pc at the distance of each IRDC. 



tifying clumps. The information in the third dimension enables 
us to disentangle clumps that are sufficiently well separated in 
velocity. For those sources with SCUBA 850/mi data, we cross- 
identify the clumps within half the NH3 beam (~ 20 ") at the 
higher SCUBA resolution (~ 14 ") and estimate the source 
size by a 2D Gaussian fit to the clumps in the SCUBA data. In 
cases where we do not have dust emission maps, we use the 2D 



Gaussian fit routine in GRAPHIC. It searches for the bright- 
est pixel in each velocity channel of the channel map across 
the NH3 (1,1) main component and fits a 2D Gaussian to de- 
termine source size. The output is then cross-checked over the 
different channels and the fit obtained for the brightest emis- 
sion is used. We find 1-5 clumps for each cloud. The source 
sizes after correcting for the beam size (after subtracting the 
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Fig. 2. NH3(1,1) (lower histogram in each panel) and (2,2) (upper histogram, offset from K for clarity) spectra toward the 
brightest position of the 9 observed sources. The brightest positions are referred to as NH3 PI in Table|3] 



Gaussian beam size in quadrature) are listed in Table [5] The 
dense cores within these IRDCs are thus not resolved with the 
40 " 100m beam. We find several secondary peaks offset by 
> lcr in most sources; offsets from the central position are given 
in Table Table [3] should be referred to for the nomenclature 
of the clumps identified in the NH3 maps and/or the SCUBA 
850 jum images hereafter. 



4.3. Line Profiles 

In general, NH3 lines are bri ghter than the mm rotational H2CO 
transitions observed bv lCarev et alJJl998h . As given in Tabled 
the (2,2) peak line intensity is on average 40% of that of the 
(1,1) line. Gaussian line profiles with extended wings have 
been reported for H2CO lines in most of the clouds JCarev et"all 
1998) but we do not observe any pronounced wings in NH3. 
Leurini et al. 2006 (submitted) find line profiles similar to those 
observed in H2CO in the mm transitions of CH3OH. 
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Fig. 4. SCUBA 850yum image of G28.34 with NH 3 (1,1) in- 
tegrated intensity as contours. The contour levels are (2cr, 4cr, 
6a...). 



4.4. Linewidth 

The width of a spectral line is an important parameter because 
it is a measure of the total kinetic energy of the cloud. The 
NH3(1,1) linewidths for all sources lie between 0.8 - 3 km s _1 . 
The H2CO linewidths are significantly larger (factor 2) than 
the NH3 linewidths in all cases except G8 1.50+0. 14, where the 
NH3 linewidths are twice as large. 

The NH3 linewidths for our sample ar e higher than those 
of NH3 cores reported in Ijiiina et alJ (H"999). which are mostly 
low mass cores. The large linewidths might be explained as 
due to clumping. Clumps with smaller linewidth but with a 
higher clump-to-clump velocity dispersion may add up to the 
observed linewidths. The sources closest to us (G79 IRDCs) 
have the smallest linewidth and belong to the Cygnus-X re- 
gion. In order to test whether the larger linewidth we observe 
at larger distance is a distance effect, we average the G79 IRDC 
linewidth over an area 5A oc (dt m ldcy g ) 2 where df m - is the dis- 
tance to a source farther than Cyg and with a much higher 
linewidth. On comparing the resultant scaled values, we find 
that the linewidths are still significantly larger for the sources 
at larger distances after accounting for the larger region in each 
beam for the more distant cores. 

To illustrate this point, let us take the example of the 
two extreme values of linewidths from our sample. Dark core 
G79.27+0.38 NH 3 P2 (at 1 kpc) has a linewidth at the bright- 
est position of 0.83 km s 1 while G28.34+0.06 NH 3 P2 (at 



4.8 kpc) has a linewidth of 2.65 km s -1 . Neither of the two 
cores have a MIR counterpart which might indicate any deeply 
embedded protostar influencing its immediate environments. 
Averaging all the emission in G79.27+0.38 over an area of 
~ 200 x 200", we obtain a linewidth of 1 .5 km s~ 1 , still smaller 
than 2.65 km s _1 . Therefore, the large linewidths derived ap- 
pears to be mainly due to the velocity dispersions within the 
beam. 

One of the assumptions made to derive the temperature is 
that the beam filling factor is the same for both inversion lines; 
that is, the two lines are emitted by the same volume of gas. The 
(2,2) linewidths are slightly larger than the (1,1) linewidths (cf. 
Table|2ji toward only some cores, suggesting that the same gas 
is not exactly traced by both lines in those sources. 

The linewidths exhibited by these sources are much larger 
than the thermal linewidth, Av t h, which for 7\i n ~ 2QK should 
be ~ 0.22 km s _I as per the relation 

Avth ~ a I , (2) 

where k is Boltzmann's constant and wnh, is the mass of the 
ammonia molecule. Average linewidths ~ 2 km s _I could be 
explained in terms of velocity gradients due to rotation of the 
cloud, or turbulent cloud movements. The linewidths seems 
to be especially high towards the cloud G28.34+0.06, where 
there are several unresolved NH3 clumps, identified in the 
SCUBA images. In the G79 IRDC the linewidth decreases from 
G79.34+0.33 NH 3 PI (1.22 ± 0.03 kms" 1 ) to G79.27+0.38 
NH 3 P3(0.82+0.03 km s' 1 ). As shown in Fig.|5] we find that the 
linewidth anti-correlates with optical depth for G28.34 while 
there is a possible positive correlation G79.27. 

G79.27+0.38 NH3 PI seems to be more quiescent than its 
surroundings. In the G79 cloud complex, there is considerable 
difference in linewidth (~ 0.5 kms 1 ) between the different 
clumps along the cloud from east to west. Thus, these clouds 
might be harb ouring objects at different stages of evolution 
dRedman et al.l20 03). 

4.5. Kinetic Temperature 

Since the optical depth is known, we derive the excitation tem- 
perature of the NH3 (1,1) inversion transition (Table [3} via the 
relation 

?ex = 7— —7 + 2-7 K, (3) 
1 - e T 

where Tmb and t represents the temperature and the optical 
depth of the of the (1,1) line. By fitting the main and the hyper- 
fine components of the (1,1) line and the main component of 
the (2,2) line, w e obtain the rotational temperature. An analyt- 
ical expression jTafalla et al.ll2004T) has been used to estimate 
the kinetic temperature from the rotational temperature. The 
kinetic temperature is given by the expression 



l-^ln[l + l.lexp(-16/T rot )]' 
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Fig. 5. Correlation plot between the NH3(1,1) optical depth and 
the NH 3 (1,1) linewidth. Bottom panel: G79.27+0.38 NH 3 PI, 
Top panel: G28.34+0.06 NH 3 P2. 

The typical temperatures for the cores range from 1 1 - 17 K 
as displayed in Table|3 In some cases, we find that the fit to the 
NH 3 (1,1) line slightly underestimates the peak intensity. 

We find that the kinetic temperature is significantly higher 
than the excitation temperature, as given in Tabled The beam 
filling factor r] is a measure of the fraction of the beam filled by 
the observed source. Assuming that the cores are in local ther- 
modynamic equlibrium (LTE), we may estimate this fraction 
as J] — Ttx/Tkin where T^n is the kinetic temperature of the 
gas. We derive filling factors of ~ 0.3 - 0.5 for all the clumps. 
These low filling factors suggest either sub-thermal excitation 
(non LTE conditions) or clumping within the beam. From our 
estimation of the sizes (Table [3}, we find most of the cores are 
unresolved with the 40" beam. Thus clumping is more likely 
to explain the small filling factors. 

We find hints of temperature gradients (inside-out) within 
the cores in three sources and a reverse gradient in one source. 
However the large error bars associated with the rotational tem- 
perature does not allow us to make a convincing case. The tem- 
perature structure will be analysed in detail in a future paper 
with observations at high angular resolution (Pillai et al. 2006 
in prep). 

There is also a variation in temperature within the different 
cores of the same cloud. In Gl 1.1 1-0.12, with a projected ex- 
tent of a few pc, the average (over the core size) temperatures 
of the southern cores are 2 - 3 K (10.9 + 1 K) lower than the 
northern cores (13.4 + 1 K). 



In G28.34+0.06, the gas temperature derived towards NH 3 
P3, 12.7 ± 1 K and NH 3 P5, 13.6 + 1 K is significantly lower 
than in the rest of the cloud. CH 3 OH observations on the PI 
position (Leurini et al. 2006 submitted) reveal a cold and a hot 
component, the latter with a small filling factor. We obtain ki- 
netic temperatures of ~ 16 K towards the peaks of NH 3 PI and 
NH 3 P2. Higher angular resolution observations are needed to 
confirm a hot component of small extent also in ammonia. 

In G19.30+0.07, the temperature at the position of peak 
NH 3 emission PI is 17.9 + 1.6 K, much warmer than the rest of 
the cloud (13.1 ± 1.0 K). 

In the G79 complex, we find that the temperature in 
G79.34+0.33 NH 3 PI is higher (14.3 ± 0.8 K) than the other 
2 cores G79.27+0.38 NH 3 PI (12.2 ± 0.8 K) & G79.27+0.38 
NH 3 P2 ( 1 1 . 1 ± 1 .2 K) in the western part of the filament. 

4.6. Velocity Structure 

The average velocity-position diagram along an axis with posi- 
tion angle of +40 ° (to align the position axis roughly with the 
direction of the filament on the sky) for the cloud Gil. 11 -0.12 
is shown in Fig. [6] At every position along the position axis, 
the NH 3 (1,1) spectra were averaged along a line perpendicu- 
lar to the position axis. There is a clear trend for the velocities 
to decrease from the south to the PI position and increase to 
the north of the NE filament. This is also seen in the channel 
maps shown in Fig. where the different clumps in the north 
and south appear at distinctly separate velocities. Given that 
the NH 3 (1,1) linewidth at all positions is < 2.5 kms -1 , this 
velocity shift (> 5 kms -1 ) between the peaks is certainly sig- 
nificant. The unresolved clump Gl 1.1 1-0.12 NH 3 P4 is more 
associated with the G 1 1 . 1 1 -0. 1 2 NH 3 P2 clump while P3 fades 
away at those velocities. 

Recently, this filament had been studied in absorption 
against the diffuse 8 pm Galactic background and in emis- 
sion from cold dust at 850 pm Jjohnstone et aljE003h . They 
model the 850 pm emission by fitting a non-magnetic isother- 
mal cylinder profile to the radial structure of the entire filament 
(south and north segment). The velocity structure observed in 
NH 3 (1,1), however, shows several distinct clumps. This sug- 
gests that the density structure of the filament cannot be rigor- 
ously described by a simple continuous cylinder. 

The channel maps for G28.34+0.06 are shown in Fig. [8] 
The two extensions of the cloud around PI and P2 differ in 
velocity by about 1.5 kms" 1 . This is the difference in veloc- 
ity for the two peaks at which they are brightest. From the 
channel maps it appears that there is a bridge between the two 
parts of the cloud connecting PI and P2. At the velocity of 
80.16 kms 1 , the three unresolved clumps, which are identi- 
fied in dust continuum emission but not clearly identified in the 
NH 3 (1,1) integrated intensity emission, are revealed. 

G33.71-0.01 shows another interesting case of distinct ve- 
locity variations across a cloud. The average velocity position 
map along the declination axis is displayed in Fig. [9] The veloc- 
ity increases from the south at ~ 103 km s _1 towards the north 
to ~ 106 kms 1 while showing a very wide weak component 
at ~ 105 kms -1 at the (0,0) position. Its detection at (0,0) is 
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Fig. 7. Channel maps for Gl 1.1 1-0.12 over the main component for the NH3(1,1) transition. 



only at the 2<x level but is also observed in NH3 (2,2). A possi- 
ble interaction with the nearby SNR G33. 6+0.1 is discussed in 
§5.2. 

4. 7. Column Density, NH 3 Abundance and Masses 

We have estimated the column densities and total masses for 
the bright, compact sources for which the SCUBA 850 fim 
maps are available, after smoothing the data to the resolution 
of 100m beam (40"). Assuming a dust opacity of K m of 1.85 
cm 2 /g at 850 fim for grains with t hick ice mantles a nd gas den- 
sity n(H) — 10 cm (Ossenkonf & Henning 1 1994 . the effec- 
tive H2 column density is (Launhardt 1996) 



JV(H 2 ) 



6.2 x 10 16 5^ eam /l 3 e 



1 .44x1 4 



zo 



(5) 



where S^ e ' dm is the flux density in Jy/beam, A the wavelength 
in fim, Z/Zo is the metallicity relative to the solar metallic- 
ity (we assume Z/Zo = 1) and 9 is the FWHM of the 100 m 
beam in arcseconds. Here we assume that the dust temperature 
Td is approximately equal to the gas temperature. We derive 
an effective H2 column density of the order of 6 - 20 x 10 22 
cm -2 toward the peaks of the dust emission. Together with the 
NH3 column density, this is used to determine the ammonia 
abundance (^nh 3 = N(NH3)/N(H2) ). We determine the peak 
effective H column density from the SCUBA map and the NH3 
column density for the same position for each clump and get 
an upper limit for the fractional NH3 abundance as listed in 
Table [3] ■ The average abundance is ~ 4 x 10~ 8 , however has a 
low value towards G79.27+0.38 NH 3 P2 (7 x 10 -9 , as listed in 
Table|3}. 
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Fig. 8. channel map for G28.34+0.06 over the main component for the NH3(1,1) transition. 
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Fig. 9. Left panel: The average velocity-position map of NH3 (1,1) along the declination axis for the cloud G33.71-0.01. Right 
panel: The spectra obtained at the corresponding position in the map (indicated by the arrows). Note the large linewidth compo- 
nent toward the centre part of the cloud and the change in LSR velocity from south to north. 
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(km/s) 



Fig. 6. The average velocity -position map of NH3 (1,1) for the 
cloud Gil. 11 -0.1 2 at a position angle of ~ 40 °. Note that 
the Y axis is the offset from the reference position, given in 
Table[T] offsets run from south-west (negative offsets) to north- 
east (positive offsets). 



5. Discussion 

5.1. Gas Temperature and Distribution 

In Gil. 11 -0.1 2, the position NH 3 P2 is a faint MSX source 
with no counterpart in the visible or N IR; it is very likel y that 
this is a heavily embedded protostar (Pill ai et alJl20 05). The 
southern clumps are colder by ~ 2 - 3 K than the northern 
clumps. Thus, these southern filaments may be at an earlier 
stage of evolution than the cores belonging to the northern fila- 
ment, where star formation activity might heat up the gas. The 
core to core variation in temperature within an IRDC is ob- 
served toward several other sources (see ^4.51 . Therefore, these 
dense cores in IRDCs are ideal sites for investigating the initial 
conditions in forming a massive star. 

The gas tempera tures we derive ge nerally agree with the 
dust temperatures of ICarev et alJ J2000l) based on submm ob- 
servations and an assumed dust emissivity index ji — 1 .75. r gas 
might be lower limits for those sources that show up as bright 
and compact objects in the submm dust emission, indicating 
the presence of a heavily embedded object, but reasonable for 
the cold gas envelope. 



The mass of the each cloud can be determined from the dust 
continuum and from NH3 emissi on. The gas mass is deriv ed 
from the dust using the relation! Launhardt & Henning 1997) 

1.44X10 4 

M(H 2 )[M ] = 1.8xlO- 14 Sy/ld " TJ ' gg[M B ] (6) 

K m (A) Z 

where S v is the integrated flux density in Jy and the other pa- 
rameters are the same as in eq.(|3. 

Our estimates for the mass and H2 column density are 
in ge neral agreement with the previous work of ICarev et alJ 
(2000) for the submillimeter cores. 

The total gas mass can be derived from the NH3 column 
density maps assuming a uniform fractional abundance of the 
molecule, for those sources without dust continuum data. The 
gas mass derived from SCUBA observations can be directly 
compared with the virial mass estimate. The virial parameter 
( Bert oldi & McKeel 19921) for a clump is defined as 



where <x is the three dimensional root mean square (r.m.s) ve- 
locity dispersion and R is the radius of the clump, and M is 
the gas mass. Note that cr = V3/81n2 x FWHM. The virial 
mass is defined as M v h = 5<x 2 R/G. For the clumps to be sta- 
ble against collapse, a ~ 1. We find an average value of 2.1 
for the virial parameter toward the individual cores. As the un- 
certainty in the dust opacity alone is a factor of 2, a ~ 1, and 
most of these cores appear to be virialised. Therefore, the core 
structure is consistent with them being supported by turbulent 
pressure without any evidence of external bounding pressures. 
G28. 34+0.06 NH 3 P4 is an extreme case, where a = 6.53. 



5.2. Physical Properties, Chemistry & Lifetime 

The excellent correlation between the NH3 line and submm 
continuum emission is consistent with the predictions of chem- 
ical models that NH3 is relatively more abunda nt in high 
density region than other molecules (Bergin 2003). The col- 
umn densities translate to extremely high extinction values of 
55 - 450 mag. The discre pancie s in NH3 and H2CO linewidths 
reported by ICarev et alJ ( 119981) (see ^4.4> indicate the differ- 
ences in the gas volume traced by the two molecules. The line 
wings seen in H2CO and CH3OH (Leurini et al. 2006 submit- 
ted) gas toward some sources is most likely high-velocity out- 
flow. The critical density of the H2CO lines studied (10 6 cm 4 ) 
is much higher than that of the NH3 inversion lines (10 4 cm 3 ). 
Hence, while physical parameters derived from NH3 are rep- 
resentative of the general, cool IRDC material, H2CO probes 
the dense gas, which might be influenced by embedded pro- 
tostars. The derived NH 3 abundance of 0.7 - 10.1 x 10~ 8 to- 
gether with the centrally condensed NH3 emission is consistent 
with the chemical mod el predicted for pre-protostellar cores 
dBergin & LangeJl997l) . These values imply that NH3 is over- 
abundant by factors of 5 - 10 relative to "normal"(= lower 
density) and less turbulent dark clou ds. In contrast, H ?CO is 
under-abundant by a factor of ~ 50 dCarev et alJl 998). Hence 
the IRDCs appear to exhibit a complex chemistry. 

We derive the flow crossing time, for the cloud to disperse 
due to its own internal motions. Perturbations within the cloud 
would disintegrate the cloud unless the propagation wavelength 
is much smaller than the distance across the cloud which it tra- 
verses. The estimated dynamical timescales of ~ 0.5 - 2 Myrs, 
provide an upper limit to the life time of these clouds. 

The virial masses for the different clumps span a wide range 
from 50 to 7880 M . As typical for intermediate to high mass 
star-forming regions, most of the masses are skewed to values 
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> 500 M Q . The virial parameter 
which are stable against collapse. 



1 for most of the clumps 



5.3. IRDCs and Other Star-forming Regions 

In a comparison of our NH3 resul ts with the potential low mass 
counterpart B68 llLai et alJ l2003). we find that the linewidths 
and the virial masses are much higher towards our cores. The 
sizes and the column density and most importantly the masses 
we derive are also significantly higher, thus making IRDCs 
potential candidates as sites of high-mass star formation. The 
(1,1) and (2,2) linewidths for IRDCS are slightly lower but the 
NH3 column densities are comparable to the m assive de nse 
cloud identified in NH 3 in NGC6334 I(N) by iKuiper et al] 
i ll 9951) . In terms of masses a nd sizes and temp eratures, our 
sources are similar to those of iGarav et al.l d2004l) . Star forma- 
tion probably has already started in some of them. 

In Table |4] we present the mean properties of the dense 
cores of our sample with those of the co res in Taurus, Per seus 
and the Orion A complex as given in Ladd et al. (1994J). It 
is evident that IRDC cores are on average highly turbulent, 
larger with much higher masses than the other cores while hav- 
ing similar temperature as the other cores. The IRDC masses 
deri ved here are for an NH3 abundance of 10~ 7 as used by 
lLadd etalJ dl994l) . But this is almost an order of magnitude 
higher than what we derive and hence the mean mass quoted in 
Table |4] will be an order of magnitude higher, several 1000 M 
instead of several 100. Compared to local dark clouds, IRDCs 
pile up significantly large amounts of mass and have supersonic 
internal motions. But how much of this mass goes into forming 
stars of low, intermediate or high mass is yet to be answered. 

iTanl J2005I) find that a group of local IRDCs have masses of 
a few 10 3 - 10 4 M and mean surface density of £ ~ 0.1 gcirT 2 . 
They find that this is 3 times the mean surface density of a Giant 
Molecular Cloud (GMC) and very similar to the values found 
in more evolved systems like the Orion Nebula Cluster (ONC). 
Subsequently, they suggest that IRDCs forming from GMCs 
are the initial conditions for star clusters. 

120051) define £ as E = M/(nR 2 ), and if we compute the 
mean surface density for our sample and the local dark clouds 
(the mean Mass and the size as given in Table @J we arrive at 
values between 0.08 - 0.4. If the mean surface density were to 
be a measure of star formation efficiency, then IRDCs have a 
higher value compared to Taurus but not very much higher than 
Perseus which is understood to be an intermediate star forming 
region. 

Studies of relatively local cluster forming regions like the 
p Ophiuchus cloud, the Serpens molecular cloud and Orion 
B molecular cloud dMotte et alJll998l iTesti & Sargeiitlll998l 
llohnstone et alJ l200ll) find that the mass spectrum of their 
cores are very similar to that of the stellar IMF. If this is in- 
deed true, then the fraction of the core mass going into form- 
ing stars would be independent of mass and the stellar IMF 
woul d mainly be determined by the cloud fragmentation pro- 
cess dBlitz & Williams! 1 19991) . If we assume a star formation 
efficiency of 30% in an IRDC core of mean mass ~ 500 M and 
adopt the standard IMF (with power law indices 2.3 + 0.3, 



1.3 + 0.5 and 0.3 ± 0.7 for masses > 0.5, 0.08 - 0.5 and 
< 0.05 M respectively), then ~ 116 stars could form in 
the core. Out of this, ~ 2 stars could be of high-mass stars 
(> 8 M Q ) while 64 would be low-mass/intermediate stars 
(0.5 < M/M < 8) while the rest would be very low mass and 
sub-stellar objects. In deed there is growing evidence of star for- 
mation in these co res dRathborne et al 120051 IPillai et all2 005. 
lOrmel et alfeool. 



5.4. IRDCs Formed by Interaction with SNRs ? 

From the velocity structure, Gil. 11-0. 12 NH3 P2 appears blue 
shifted while the clump NH3 PI towards the north and NH3 P4 
towards the south appears redshifted relative to the LSR veloc- 
ity. This is thus not a case of smooth velocity gradient along 
the filament. It might be possible to explain the observed ve- 
locity structure, if we assume that the entire filament seen in 
projection is part of a unbound system where PI and P4 lie at 
the same distance along the line-of-sight (l.o.s) at diametrically 
opposite ends while P2 is further in the front. A massive wind- 
driven process which might have taken place in the close vicin- 
ity of the cloud could explain such a structure. Recent wide 
field MIR images of this region released from Spitzer show that 
this cloud has a filamentary concave structure spread over sev- 
eral par^e£s_wifiia^significant density enhancement in the cen- 
ter dMenten. Pillai. & Wvrow ski 2005). The morphology bears 
remarkable similarity to the structures predicted by 3-D numer- 
ical calculations simulating the impact of a planar shock front 
on an isolated globule (Bossl ll995l) . The best known case of 
a ongoing SN-cloud interaction is in IC 443, where a 10 4 yr 
old SNR in the GEM OBI association with a shock speed of 
40 km s _I impacts the cloud resulting in highly excited molec- 
ular gas with very broad linewidths. However, a weaker shock 
(likely from a more evolved SNR) would result in temporary 
distortion and compression of the molecular cloud, followed 
by rebound to a equilibrium. An evolved SNR at an age of 
10 5 yr with a shock speed o f 100 km s - 1 , is one of the three 
likely weak shock waves that Boss ( 1995) propose to influence 
a cloud without destroyi ng it. Such a s hock wave would have 
already traversed 25 pc.l Brogan et alJ d2004h recently discov- 
ered the supernova remnant Gil .03-0.05 which is within 7' of 
the cloud. They claim that the SNR is not young. 

The SNR has a shell-like appearance with a diameter of 
~ 8 ' and very w eak emission at centimeter wavelengths. 
iBro gan et al] J2004 find that a pulsar PSRJ1809-1917 is lo- 
cated about 8.5' from the SNR and is probably associated 
with the SNR. Considering the respective uncertainties of their 
distance estimates, the pulsar (4+1 kpc) and Gil. 11 -0.1 2 
(3.6 + 0.5 kpc) appear at or near the same distance. 

In G33.71-0.01, we find another where the cloud seems to 
have undergone a shock at the centre (see Fig. Note that 
there is a large linewidth towards the central core of this fil- 
ament. The (0,0) velocity is also different from the cores to- 
ward the south and north segment of the filament. The SNR 
G3 3.6+0. 1 is w ithin 5' of this core and previous observations 
by iGreenl dl989l) shows an unusually broad OH absorption 
feature near 105 kms 1 towards this SNR. The HCO + and 
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12 CO J = 1 - observations bv lGreen & Dewdnevl dl992l) re- 
veal material shocked by the interaction of the SNR with the 
adjacent molecular cloud. The 2cr detection of a large linewidth 
feature Av > 8 kms -1 centered at 105 kms -1 in NH3 (1,1) and 
(2,2) might also be a sign of interaction of the cloud with the 
SNR. 

All IRDCs except those in the Cygnus-X region have a 
SNR in their vicinity, although an association, which would 
require observations with shock tracers, cannot be verified yet. 
To estimate the chance occurance of an IRDC and SNR o n the 
sky, we sample ~ 17.6 s quare degrees in which fcireerJl2002h 
and lBrogan et alJ (E004) have identified a total of 22 remnants. 
The likelihood of a chance positional occurance within 7' of 
an IRDC is small (< 0.05) given a SNR surface density of 
1.25 remnants per square degree. It is very likely that the SNRs 
are physically associated with the IRDCs Gil. 11 and G33.71- 
0.01. 



5.5. IRDCs in the Framework of an Evolutionary 
Sequence 

To obtain general statistics on IRDCs we compare in Fig. EH 
the core gas properties of our sample o f 9 sources with the 
cores presented in the Ijiiina et all i 19991) NH3 database. The 
database consists of 264 dense cores, with and without asso- 
ciations with young stellar objects. For the large linewidths 
of IRDCs (Av > 0.8 kms 1 ), there is a distinct trend for the 
IRDCs to be colder relative to the cores from the database. The 
NH3 column density is also higher for the IRDCs, surpassed 
only by a few high column densities from low mass cores ob- 
served with the high angular resolution of the VLA, which ev- 
idently sees more of the core interiors that are possibly associ- 
ated with YSOs. 

In Fig. II llwe present a simil ar comparison with a sample of 
UCHn regions from the Wood & Churchwell ( 1989) catalogue 
and the sources from the iBeutheretall (2002) study of high- 
mass protostellar objects (HMPOs). While recent stud ies reveal 
that HMPOs are in a pre- UCHn region phase ( Mol inari et alJ 
I2002L lBeiitheretalJ |2002). the nature of the stage preceding 
HMPOs has not been studied in detail. IRDCs are ideal candi- 
dates for this pre-HMPO stage. In Fig.^2 we present proper- 
ties of source samples believed to cover the earliest phases of 
massive star formation based on NH3 observations. There is no 
clear trend in linewidths, however IRDCs have a significantly 
lower average linewidth (Av = 1.51 kms _I ) than HMPOs 
(Av = 2.05 kms" 1 ) and UCHn regions (Av = 2.52 kms" 1 ). 
There is a distinct temperature trend from the low tempera- 
tures of the IRDCs (T = 13.9 K) to increasing temperatures 
for the IRAS selected high-mass objects (T = 18.34 K) and 
high temperatures for the objects associated with UCHn re- 
gions (T = 22.6 K). Temperature and linewidths must be un- 
derstood as averages over the core and parts of the envelope 
since the beam is 40". 

There is no clear trend in column density. However, on av- 
erage the column densities of IRDCs are high compared to the 
other two samples. We interpret the clear trend in temperature 
and the tentative differences in NH3 linewidths and column 
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Fig. 10. Bottom panel: The correlation plot of kinetic tempera- 
ture with the linewidths. The black filled squares indicate the 
IRDCs and the grey unfilled squares the cores from the NH3 
database. IRDCs have larger velocity dispersions while the ro- 
tational temperatures are comparable. Top panel:Tbe correla- 
tion plot of NH3 column density with the linewidths shows that 
the average column density of IRDCs are higher. 



densities (N[NH3]) as possible manifestation of an evolution- 
ary sequence. Starless cores on the verge of star formation are 
expected to be cold (T<20 K), to have high column densities 
and smaller linewidths. The temperature and the linewidth in 
a core will increase after the formation of an embedded proto- 
star (via radiative heating and injection of turbulence through 
outflows and winds), while the envelope will be dispersed via 
outflows and winds and thus column densities will decrease 
with time. 

One would expect the highest ammonia column densities 
for the hot sources due to evaporation of ammonia from the 
grains although in some cases IRDCs have higher column den- 
sities. However, this effect might be very localised to the hot 
cores (0. 1 pc and smaller) and not be true for the larger scale 
emission. Additionally some molecular material might already 
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Fig. 11. Bottom panel.The correlation plot of kinetic tem- 
perature with the linewidth. The black filled squares indi- 
cate the IRDCs,the crosses the sources from the Woo d and 
Churc hwell catalogue and the unfilled squares iBeuther et alJ 
sources. IRDCs are colder and on average have lower 
linewidths. 7b/? panel:The correlation plot of NH3 column den- 
sity with the linewidths. 

be dispersed in the hot sources by the interaction of the young 
OB clusters with their environment. 



6. Conclusions 

In this paper, we discussed ammonia observations of a selected 
set of IRDCs and the derived physical properties. Our results 
are as summarised below. 

The ammonia emission correlates very well with MIR ab- 
sorption and ammonia peaks distinctly coincide with dust con- 
tinuum peaks. Several cores are detected within the clouds with 
deconvolved sizes smaller than the 40" FWHM beam size. We 
can constrain the average gas temperature to between 10 and 
20 K. 

We observe high linewidths (1 < Av/km s _1 < 3.5), hence 
turbulence plays an important role in the stability of an IRDC. 
There are significant velocity gradients observed between the 



cores. The effect of external shock/outflow tracers , on the gas 
kinematics is suggestive in some cases, but needs to be investi- 
gated further. 

The column densities translate to extremely high Ay values 
(55 - 450 mag), therefore any active star formation would be 
heavily embedded. The total cloud gas mass derived from the 
NH3 data ranges from 10 3 - 10 4 M . The virial parameter is 
~ 1 for most of the clumps, and the cores appear to be stable 
against gravitational collapse. As a result IRDCs are potential 
sites for star formation. If we were to adopt the stellar IMF and 
a star formation efficiency of 30%, then every IRDC core could 
fragment to form > 100 stars, with at least two high mass stars 
(> 8 M ). 

The fractional abundance of NH3 (relative to H2) is 0.7 - 
10. Ox 10~ 8 . This, together with the excellent correlation in mor- 
phology of the dust and gas, is consistent with the ti me de pen- 
dent chemical model for NH3 of Bergin & Langer ( 1997) and 
implies that NH3 remains undepleted. The derived abundance 
is a factor 5-10 larger than that observed in local dark clouds 
while H2CO is underabundant by a factor of ~ 50. Hence, the 
chemistry governing these IRDCs might be complex and could 
be different from other parts of the dense ISM. 

The time scales we derive for the clouds to disperse due 
to their own internal motions of a few Myrs provides an up- 
per limit to the life time of these clouds. We suggest that SNRs 
might be the trigerring mechanism responsible for the forma- 
tion of an IRDC. 

The comparison of the physical properties from ammonia 
of our IRDCs sample with other source samples - HMPO's 
and UCHns strongly suggests that most of these IRDCs are the 
most likely candidates for pre-protostellar cores of massive star 
formation. 

Acknowledgements. This research has made use of the NASA/ IPAC 
Infrared Science Archive, which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with the 
National Aeronautics and Space Administration. S.J.Carey acknowl- 
edges support from a NASA Long Term Space Astrophysics grant. 
T.Pillai was supported for this research through a stipend from the 
International Max Planck Research School (IMPRS) for Radio and 
Infrared Astronomy at the University of Bonn. T.Pillai thanks her 
colleague J.Kauffmann for useful discussions and comments on the 
manuscript. 

References 

Bachiller, R., Guilloteau, S., & Kahane, C. 1987, A&A, 173, 324 

Bergin, E. A. 2003, in SFChem 2002: Chemistry as a Diagnostic of 
Star Formation, proceedings of a conference held August 21-23, 
2002 at University of Waterloo, Waterloo, Ontario, Canada N2L 
3G1. Edited by Charles L. Curry and Michel Fich. NRC Press, 
Ottawa, Canada, 2003, p. 63., 63-+ 

Bergin, E. A. & Langer, W. D. 1997, ApJ, 486, 316 

Bertoldi, F. & McKee, C. F. 1992, ApJ, 395, 140 

Beuther, H., Schilke, P., Menten, K. M., et al. 2002, ApJ, 566, 945 

Blitz, L. & Williams, J. P. 1999, in NATO ASIC Proc. 540: The Origin 
of Stars and Planetary Systems, 3 — v 

Boss, A. P. 1995, ApJ, 439, 224 

Brogan, C. L., Devine, K. E., Lazio, T. J., et al. 2004, AJ, 127, 355 
Carey, S. J., Clark, F. O., Egan, M. P., et al. 1998, ApJ, 508, 721 



14 T. Pillai, F. Wyrowski, S.J. Carey and K.M. Menten: Ammonia in Infrared Dark Clouds 

Table 1. List of IRDCs Observed in NH 3 (1,1) and (2,2) 



Name 


R.A.(J2000) 


Dec.(J2000) 


Vlsr [kms '] 


D [kpc] 


G11.11-0.12P1 


18:10:29.27 


-19:22:40.3 


29.2 


3.6 


G19.30+0.07 


18:25:56.78 


-12:04:25.0 


26.3 


2.2 


G24.72-0.75 


18:36:21.07 


-07:41:37.7 


56.4 


1 


G24.63+0.15 


18:35:40.44 


-07:18:42.3 


54.2 


3.6 


G28.34+0.06 PI 


18:42:50.9 


-04:03:14 


78.4 


4.8 


G28.34+0.06 P2 


18:42:52.4 


-03:59:54 


78.4 


4.8 


G33.71-0.01 


18:52:53.81 


+00:41:06.4 


104.2 


7.2 


G79.27+0.38 


20:31:59.61 


+40:18:26.4 


1.2 


1 


G79.34+0.33 


20:32:21.803 


+40:20:08.00 


0.1 


1 


G81.50+0.14 


20:40:08.29 


+41:56:26.4 


8.7 


1.3 



Notes: Columns are name, right ascension, declination, LSR velocity, distance. Positions and Distances are taken from lCarev etai]ll998h . The 
coordinates correspond to the reference positions of the maps. 
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Table 2. NH 3 (1,1) and (2,2) Map Results: Peak Position. 



Source 


Transition 


Vlsr 


Tmb 


FWHM 


^"main 


1 a- 








[km s~'] 


[K] 


[kms -1 ] 




[Kkms" 1 ] 


Gil. 11-0.12 NH 3 PI 


1- 


-1 


30.44 (±0.02) 


5.47 (±0.96) 


1.27 (±0.06) 


3.52 (±0.39) 


0.42 




2- 


-2 




1.7/ v^znu. J^} 








G19.30+0.07 NH 3 PI 


1- 


-1 


26.45 (±0.03) 


5.17 (±0.8) 


1.88 (±0.08) 


1.88 (±0.26) 


0.53 




2- 


-2 


26.64 (±0.06) 


2.86 (±0.27) 


1.57 (±0.17) 






G24.72-0.75 NH 3 PI 


1- 


-1 


58.37 (±0.07) 


2.76 (±0.45) 


2.57 (±0.12) 


2.81 (±0.47) 


0.25 




2- 


-2 


58.47 (±0.12) 


1.99 (±0.23) 


2.97 (±0.24) 






G24.63+0.15 


1- 


-1 


53.11 (±0.04) 


2.40 (±0.34) 


1.7 (±0.09) 


2.36 (±0.43) 


0.16 




2- 


-2 


52.89 (±0.09) 


0.98 (±0.12) 


2.42 (±0.24) 






G28.34+0.06 NH 3 PI 


1- 


-1 


79.14 (±0.02) 


2.88 (±0.42) 


2.67 (±0.04) 


2.00(0.11) 


0.16 




2- 


-2 


79.15 (±0.04) 


1.37 (±0.05) 


3.29 (±0.11) 






G33.71-0.01 NH 3 PI 


1- 


-1 


104.31 (±0.09) 


2.52 (±0.43) 


2.64 (±0.2) 


1.5 (±0.52) 


0.12 




2- 


-2 


104.4 (±0.17) 


1.14 (±0.22) 


2.40 (±0.35) 






G79.27+0.38 NH 3 PI 


1- 


-1 


1.65 (±0.02) 


2.28 (±0.38) 


1.04 (±0.05) 


2.49 (±0.27) 


0.31 




2- 


-2 


1.63 (±0.05) 


0.61 (±0.08) 


0.98 (±0.14) 






G79.34+0.33 NH 3 PI 


1- 


-1 


0.31 (±0.01) 


2.53 (±0.30) 


1.23 (±0.03) 


1.45 (±0.13) 


0.19 




2- 


-2 


0.25 (±0.04) 


0.79 (±0.07) 


1.29 (±0.13) 






G81.50+0.14 


1- 


-1 


-4.53 (±0.05) 


2.21 (±0.35) 


1.35 (±0.15) 


0.37 (±0.46) 


0.27 




2- 


-2 


-4.69 (±0.12) 


0.67 (±0.21) 


1.00 (±0.34) 







Notes: Columns are name, NH 3 (J,K) transition, LSR velocity, (1,1) main beam brightness temperature, full linewidth at half maximum, main group optical depth and the 1 sigma noise level in 

the NH 3 (1,1) integrated intensity map Q. The error estimates are given in brackets. 



T3 



| 

8 
s 



o 



> 

3 
3 
5 

3 



o 



Q 

5* 
c 

p. 



Table 3. Physical Properties of Observed IRDCs. 



Source 


Offsets 


Tex 


Trot 


Tktn 


A<(NH 3 ) 


N(H + H 2 ) 


*NH 3 


Mass 


Size 


M vir 


a 


Mnh 3 




[arcsecs] 


[K] 


[K] 


[K] 


[10 14 crrr2] 


[10 22 crrr2] 


[io- s ] 


Me 


[arcsecs] 


[M ] 




[10 4 M G ] 


G11.11-0.12NH, PI 


(80,40) 


7.6 (±0.7) 


12.3 (±1.2) 


12.7 


41.5 (± 6.4) 


8.4 


4.9 


485 


27 


468 


1.0 


18.1 


Gil. 11-0.12 NH 3 P2 


(0,0) 


6.6 (±0.4) 


13.4 (±1.0) 


13.8 


30.5 (± 3.3) 


10.2 


2.9 


172 


13 


268 


1.6 




Gil. 11-0.12 NH 3 P3 


(-300,-325) 


7.1 (±0.3) 


14.1 (±9.7) 


14.5 


60.4 (± 11.4) 


8.7 


6.9 


647 


32 


860 


1.3 




Gil. 11-0.12 NH 3 P4 


(-300,-400) 


5.7 (±0.5) 


10.5 (±1.0) 


10.8 


50.9 (± 10.6) 


16.1 


3.2 


788 


25 


555 


0.7 




G19.30+0.07 NH 3 PI" 


(20,20) 


8.2 (±0.9) 


17.6 (±2.0) 


18.4 


25.9 (± 3.2) 








38 


893 




0.9 


G19.30+0.07NH 3 P2 " 


(-60,-20) 


5.9 (±0.7) 


13.8 (±1.8) 


14.3 


33.3 (± 6.9) 








45 


823 






G24.72-0.17NH 3 PI " 


(150,220) 


6.6 (±0.6) 


17.4 (±2.2) 


18.2 


44.3 (± 5.7) 








34 


2450 




1.4 


G24.63+0.15" 


(0,0) 


5.0 (±0.5) 


14.3 (±1.4) 


14.8 


24.1 (± 4.4) 








43 


1358 




0.7 


G28. 34+0.06 NH 3 PI 


(0,0) 


5.7 (±0.2) 


16.0 (±1.2) 


16.6 


32.1 (± 2.2) 


5.8 


5.5 


904 


22 


2280 


2.5 


25.4 


G28. 34+0.06 NH 3 P2 


(20,200) 


5.5 (±0.2) 


15.4 (±1.4) 


16.0 


38.6 (± 3.2) 


32.7 


1.2 


2310 


45 


2928 


1.3 




G28.34+0.06 NH, P3 


(-60,-60) 


4.8 (±0.2) 


12.8 (±1.2) 


13.2 


47.2 (± 6.1) 


5.9 


8 


374 


21 


1398 


3.7 




G28. 34+0.06 NH 3 P4 


(-30,50) 


4.8 (±0.2) 


15.9 (±2.2) 


16.4 


37.4 (± 4.9) 


3.7 


10.1 


147 


16 


960 


6.5 




G28.34+0.06 NH 3 P5 


(50,40) 


5.6 (±0.3) 


13.6 (±1.3) 


14.1 


29.6 (± 3.8) 


6.3 


4.7 


528 


24 


1378 


2.6 




G33.71-0.01 NH 3 PI " 


(-40,-200) 


5.6 (±1.1) 


16.6 (±2.3) 


17.2 


23.0 (± 6.6) 








52 


7880 




15.6 


G33.71-0.01 NH 3 P2 " 


(60,100) 


5.3 (±0.3) 


17.1 (±2.0) 


17.8 


29.0 (± 3.3) 








60 


4928 






G79.27+0.38 NH 3 PI 


(-20,0) 


4.8 (±0.2) 


12.4 (±0.9) 


12.8 


20.5 (± 2.5) 


14.2 


1.4 


84 


32 


155 


1.8 


0.4 


G79.27+0.38 NH 3 P2 


(-140,40) 


6.0 (±0.7) 


11.2 (±1.4) 


11.5 


11.3 (± 2.7) 


12.1 


0.9 


55 


25 


53 


1.0 




G79.27+0.38 NH 3 P3 


(-250,70) 


3.3 (±2.2) 


11.7 (±1.6) 


12.0 


10.0 (± 25.0) 


6.5 


1.5 


19 


20 


40 


2.1 




G79.34+0.33 NH 3 PI 


(0,0) 


5.8 (±0.3) 


14.1 (±0.8) 


14.6 


11.9 (± 1.0) 


17.5 


0.7 


143 


39 


178 


1.2 


0.2 


G8 1.50+0. 14" 


(260,-375) 


5.6 (±1) 


14.9 (±1.3) 


15.4 


7.6 (± 2.0) 








78 


570 




0.1 



Notes: Columns are name, Offset Position from reference, excitation temperature, NH3 rotational temperature, kinetic temperature, NH3 column density, H2 column density, NH3 abundance, gas mass estimated from SCUBA data, size of the cores, 
virial masses from NH3 the alpha parameter defined as in total cloud mass from Nf^The formal errors are given in brackets. " are those sources for which the sizes have been derived from NH3 data. The H2 column densities are estimated from 850 



T. Pillai, F. Wyrowski, S.J. Carey and K.M. Menten: Ammonia in Infrared Dark Clouds 
Table 4. Mean Values of NH 3 Core Samples. 
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Sample 


FWHM(kms') 




size(pc) 


M NH , (M ) 




Taurus 


0.33 


10 


0.06 


1.1 


0.08 


Perseus 


0.55 


13 


0.12 


9 


0.16 


Orion 


1.12 


17 


0.15 


21 


0.24 


IRDC cores 


1.7 


15 


0.57 


492 


0.40 



Notes: Columns are the source sample, the mean FWHM of the NH3 (1,1) line, kinetic temperature, size of the core, mass as determined from 
NH3 assuming a fractional abun dance of NH3 rela tive to H2 and the mean surface density. The mean values for dense cores in Taurus, Perseus 
and Orion have been taken from lLadd et aljfl994h . 



